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Abstract 

By studying the Paschos-Wolfenstein (PW) ratio of deep inelastic uYe and z/Fe 
scattering cross sections, we show that it should be possible to explain the NuTeV 
sin^ anomaly with quite conventional physics, by introducing mutually different 
nuclear modifications for the valence-w and valence-c? quark distributions of the protons 
in iron. Keeping the EKS98 nuclear modifications for uy + (iy as a baseline, we find 
that some 20-30 % nuclear modifications to the uy and dy distributions account for the 
change induced in the PW ratio by the NuTeV-suggested increase Asin^6'w = 0.005. 
We show that introduction of such nuclear modifications in uy and dy individually, 
does not lead into contradiction with the present global DGLAP analyses of the nuclear 
parton distributions, where deep inelastic lA scattering data and Drell-Yan dilepton 
data from pA collisions are used as constraints. We thus suggest that the NuTeV result 
serves as an important further constraint in pinning down the nuclear effects of the 
bound nucleon PDFs. We also predict that if the NuTeV anomaly is explained by this 
mechanism, the NOMAD experiment should see an increase in the weak mixing angle 
quite close to the NuTeV result. 



^ kari . eskola@phys .jyu.fi 
^hannu. paukkunen@phys.jyu.fi 



1 Introduction 



A few years ago, the NuTeV collaboration at Fermilab reported that their measure- 
ments in deep inelastic neutrino-iron scattering indicate the value of Weinberg weak 
mixing angle sin^ = 0.2277 ± 0.0013(stat) ± 0.0009(syst) p. This result was sur- 
prising, being about 3a above the world average sin^ 6^^^ = 0.2227 ± 0.00037 A 
number of possible solutions to this deviation — the 'NuTeV anomaly' — has been 
proposed, ranging from conventional explanations within the Standard Model (SM) to 
more exotic ones requiring novel, beyond-SM, physics. For reviews, see Refs. II]. 
Today, the NuTeV anomaly is still an open problem calling for an answer. 

As intriguing it would be to see new physics appearing, it is obviously crucial to first 
investigate in detail whether the explanation could be hidden in the features of parton 
distribution functions (PDFs) of the free and bound nucleons which are not yet fully 
known. For instance, the possible local asymmetry between strange and antistrange 
quark distributions, S~ = J dxx[s{x) — s{x)] 7^ [SEllZj, and isospin violating PDFs 
[7j may both be viable candidates for explaining the anomaly. Also the PDFs of 
bound nucleons, the nuclear parton distributions (nPDFs), have been considered as 
a source for the NuTeV anomaly [S]; for a review, see 0. In particular the nuclear 
effects in the valence quark distributions uy and dy individually, have been thoroughly 
investigated in a global DGLAP analysis where the deep inelastic lA scattering (DIS) 
cross sections, the Drell-Yan (DY) dilepton cross sections in pA collisions and sum rules 
are used as constraints CHI IHl 1121 • Typically, however, the effects of the nPDFs, 
especially the differences of nuclear modifications in uy relative to those in dy, have 
been found to be too small to account for the NuTeV anomaly. 

In this paper, we consider the NuTeV anomaly and the nPDFs from a new angle, 
taking the NuTeV result as an additional data constraint for disentangling the nuclear 
effects of uy and dy distributions in nucleons (protons) of iron. First, adopting the 
EKS98 global analysis of nPDFs as a baseline, we examine in a rough but trans- 
parent way, how large nuclear effects relative to the EKS98 modification of uy + dy 
in a proton of iron would be needed to fully account for the NuTeV anomaly. Nuclear 
modifications for uy and dy deviating from the EKS98 at most on a ~ 30 % level (at 
the average NuTeV scale = 20.5 GeV^) turn out to be sufficient. 

Second, armed with the nuclear modifications for uy and dy individually, we re- 
compute the DIS and DY cross sections, the data on which constrain the EKS98 
modifications. Relative to the cross sections computed with EKS98, we find less than 
one percent changes for iron. This insensitivity shows that there can be fairly large 
differences between the nuclear effects in uy and dy which the standard global nPDF 
analyses [TSl El HI] cannot disentangle. It also demonstrates the possible role of the 
NuTeV result as a new, orthogonal, constraint in pinning down the individual nuclear 
effects in the uy and dy distributions of bound nucleons. Not forgetting the other pos- 
sible PDF-based explanations of the NuTeV anomaly, we thus suggest that a large part, 
if not all, of the NuTeV anomaly can be explained by introducing mutually different 
nuclear modifications for valence u and d quark distributions of protons in iron. 
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Third, as a consequence and a test of our hypothesis, we apply the obtained nuclear 
valence quark modifications to the neutrino cross sections measured in the NOMAD 
experiment at CERN |[14J- We predict that the nonidentical nuclear effects in uy and 
dv induce a change in the ratio of neutral-to-charged current cross sections which, if 
identical nuclear effects were used, would correspond to an increase in sin^ 9w quite 
close to that obtained by NuTeV. 

The rest of the paper is organized as follows: In Sec. 2, we define the neutrino cross 
sections and structure functions in terms of PDFs, as well as the Paschos-Wolfenstein 
(PW) ratio we study. In Sec. 3, we define the nPDFs and show how the nuclear uy 
and dv distributions enter the PW ratio. We also show analytically the direction to 
which the nuclear effects in uy and dy should deviate relative to the modification of 
uy + dy from EKS98. The results we obtain for individual nuclear modifications of uy 
and dy are presented in Sec. 4 both for the case where the NuTeV anomaly is fully 
explained and for the case where a third of the anomaly is accounted for. Section 5 
contains the discussion on the DIS and DY cross sections. The NOMAD prediction is 
made in Sec. 6, and a brief outlook is given in Sec. 7. 



2 Cross sections and Paschos-Wolfenstein ratio 

We define the framework through the following well-known charged current (CC) and 
neutral current (NC) DIS cross sections of neutrinos and antineutrinos off a nucleus. 
In particular, we shall discuss the Paschos-Wolfenstein (PW) ratio ^H]- In terms of 
double differential cross sections, the PW ratio is 

_ / dxdQ"^ — da'2^'^ / dxdQ"^ 

' da^/^'^ / dxdQ"^ — da^/"'^ / dxdQ'^'' 

where x and are the standard Bjorken DIS variables, and A denotes the target. 
In studying the consequences of the NuTeV sin^ 6w anomaly quantitatively, we shall 
focus on the PW ratio of x-integrated cross sections at a fixed scale Q^, 



daf^'/dQ' - da'/''/dQ' 
dcr'^'^'" / dQ"^ — da^jf^ /dQ"^ 



The neutral current cross section of a deep inelastic scattering of a neutrino (an- 
tineutrino) off a proton is jTHj 



dxdQ^ ~ ~\Ml + Qy 8x 



xy'F!^^'^'''^{x,Q') 

(3) 



where y is one of the standard Bjorken DIS variables and E is the (anti)neutrino beam 
energy. For these DIS variables = 2xymE, where m is the mass of a nucleon. We 
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shall use here m = 940 MeV both for protons and for neutrons. The Z-boson mass 
is Mz = 91.1876 GeV, and the Fermi coupling constant Gp = 1.16637 ■ 10"^ GeV^^ 
taken from [TB] . The structure functions are in the leading order expressed in terms of 
the quark and antiquark number densities in the free proton as 



Fr'^^(x,g2) = Ff^^(x,g2)= Y: iRl + LlMx,Q') + q{x,Q')] (4) 



q=u,d,s,c,b 



Fr"(x,Q^) = Fr"(x,Q^)=2xFr^(x,Q^) (5) 
F',^''%x, Q') = F3^'^^(x, Q')=2Y: (Ll - Rl)qyix, Q% (6) 

q=u,d 

where, as usual, Lg = — 2eqXw and Rq = —2eqXw- The weak mixing angle is 
xw = sin^ 9w, the quark's electric charge is and its third component of the weak- 
isospin is t^. All quarks are treated as massless here, as is usually the case in the 
DGLAP-evolved PDFs we shall apply below. The valence quark distribution is defined 
as qv{x, Q"^) = q{x, Q"^) - q{x, Q"^). 

The charged current deep inelastic cross section of a neutrino (antineutrino) scat- 
tering off a proton can be expressed in a similar manner [TH] . 



dxdQ"^ 









/ 2x . 



xy'F^^'^'^^{x,Q') 



For the W boson mass we take My/ = 80.425 GeV ^Hl- The charged- current structure 
functions are given by 
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Fr^^(x,Q2) = T.l^{xq.,Q') + q'{xq,Q')]\VSiu\ 

qq' 
qq' 

= 2j2[xq,q{xq,,Q') + XqY{xq,Q')]\VSiM\^ (10) 

qq' 

F2^'^^(x,g2) = 2Y.[xq'q{xq>,Q')+Xqq'{xq,Q')]\VSiu\^ (11) 

qq' 

F3^'^^(x,Q2) = 2Y.[qixq,,Q')-Y{xq,Q')]\VSL\\ (12) 

qq' 

F3^'^^(x,g2) = 2Y.[q'{xq,Q')-q{Xq,,Q')]\V^iu\\ (13) 

qq' 

where the quark flavour indices run through q = d,s,b and q' = u, c. The momentum 
fraction Xg = x{l + ^f) accounts for the mass of the produced heavy quark. We 
consider u and d quarks massless and take rric = 1.2 GeV and nib = 4.3 GeV. We 
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also set q{xq,Q'^) = when Xq > 1. In the initial state, the quark masses are ignored 
here, as is the case in the scale evolution of the PDFs, too. The Cabibbo-Kobayashi- 
Maskawa mixing matrix elements are denoted by Vq^^j. In what follows, we shall 
use ^c'km = 0.97458, I^ckm = 0-224, V^^km = 0.224, V^^km = 0.9738, ^ckm = 0.003, 
^CKM = 0.04, taken from 

For forming the Paschos-Wolfenstein ratios in Eqs. (0) and (j21), we form the differ- 
ences of the above neutrino and antineutrino cross sections. 



.(15) 



dxdQ^ dxdQ^ vr KM^ + QV 8x 2' 

xy^AF^'^ix^Q' 



dxdQ^ dxdQ^ n \M^ + Qy 2x 



where we define 



T.bvi^i'^Q")-9vi^.^Q")WcKM\^ (16) 



gg 



AF,^^{x,Q-') ^ F2^'--(x,Q^)-Fr"(x,Q^ 

= '^J2l^g'(lv{xq',Q'^) - Xqq'v{xq,Q^)]\V^iM\'^, (17) 
gg' 

SF3^^(x,g^) = F3-^^(x,Q^) + F|'^^(x,Q^) 



= 2 5: [gy(x,,, Q') + q'yixq, Q')]\VSiu\^ 
gg' 

where again the quark flavour indices run through q = d,s,b and q' = u,c. 



3 nPDFs and R 



A 



We define the quark PDFs Qa{x, Q"^) in a nucleus of mass number A and proton num- 
ber Z as in the EKS98 global DGLAP analysis Denoting the number density 
distribution of the quark flavour g in a bound proton by qP^'^(x,Q'^), and the corre- 
sponding parton distribution function in a bound neutron by g"/^(a;, Q^), the average 
number density of a quark q in the nucleus is 



Qa{x, Q'') = ZqP/^{x, Q'') + {A- Z)g"/^(x, Q 



n/Ai 



(19) 



Since we shall consider isospin effects below, it is useful to separate the symmetric and 
antisymmetric parts in the PDFs in the bound proton and neutron as 



A r 
2 



.n/A( 



+ 



A. 



g"/^(x,g^)-gP/^(x,g^) , (20) 
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where the neutron excess is denoted as = A — 2Z. In particular, as in EKS98, we 
define the nuclear PDFs to be those of the bound protons, 



g^(x,g2) = gP/^(x,Q2), (21) 

and define their nuclear modifications relative to the corresponding PDFs in the free 
proton, 

qA{x,Q^)^R^{x,Q^)q{x,Q^). (22) 

Again as in the EKS98 analysis [12], we assume that isospin symmetry between 
the bound protons and neutrons holds for arbitrary A, i.e. u^/^ = d^^^, d'^^^ = npl^^ 
and similarly for u and d. For mirror symmetric nuclei (which of course includes 
isoscalars), this gives the standard isospin symmetric way of treating the PDFs in 
protons and neutrons. Thus, in the EKS98-framework adopted here, the isospin effects 
in the average quark distributions Qa are always proportional to the neutron excess 
IS. A and to the difference between the u{u) and did) quark(antiquark) distributions in 
a bound proton, while the possible difference between e.g. the u quark distribution 
in a proton of a nucleus and the u quark distribution in a proton of its mirror 
nucleus %A is neglected. In other words, the EKS98 nuclear effects in the bound 
proton PDFs depend only on A, not on Z. The scope of the presently existing data 
does not allow for a more detailed study of these effects. In practice, however, since 
we are dealing only with close-to-isoscalar nuclei, such as iron {A = 5Q, Z = 26) here, 
the isospin-symmetry-related approximation made can be expected to be a very good 
one. 

In particular, below we shall discuss the average valence quark modification, 

^-«y(x,Q2) + rf^(^,g2)' i^^) 
in terms of individual modifications for the the valence u quarks and valence d quarks. 



Only two of these ratios are independent, so that e.g. 

^ «.(..g.)M£:g±|P« _ KA.qY^. (26) 



Analogously to the global DGLAP fits of the free proton PDFs, in the global 
DGLAP fits for the nPDFs O El CHI, the modification ratios Rf{x,Ql), i = g,q 
at an initial scale Ql are iteratively determined based on fits to the deep inelastic lA 
scattering data and the Drell-Yan dilepton data from pA collisions. Also conservation 



5 



of baryon number, charge and momentum are accounted for. The amount and pre- 
cision of the existing data do not, however, allow for a precise determination of the 
modifications R^^ and individually but for simplicity in [T31 [13 HH] it is assumed 
that at the initial scale R^^ = R^^ = Ry. Once this assumption is made, the DGLAP 
scale evolution does not cause large deviations from it, and to a good approximation 
R^^ = = R^ holds for scales > Ql as well, see EKS98 [13 • The most recent 
efforts [HI [ini [m [12] to disentangle the nuclear effects in uy and dy individually based 
on DIS and DY data do not lend support to significant differences between R^^ and 
either. 

Below, we shall study to what extent the NuTeV sin^ O-^y anomaly could be explained 
by releasing this approximation, i.e. by considering mutually different R^^ and R^^ 
but keeping their relative magnitude such that the EKS98-value of average valence 
modification, i?y, is always reproduced, as in Eq. ^I^ . In particular, we shall explore 
an extreme case where the whole NuTeV anomaly is accounted for by changes in R^ 
and Rdy, and study the possible consequences for the global DGLAP fits. 

To demonstrate the nuclear effects in the PW ratio and to see the system- 

atics expected, we first consider the ratio -R^(x, Q^) in an approximation, where only u 
and d quarks are active, V^km = virtuahties are modest, <^ M|, M^, and beam 
energy is so high that the term mxy/{2E) in Eqs. Q and ((7|) can be ignored. Using 
the nPDF definitions above, we obtain for the neutral current case 



and for the charged current case 



+ Aa^xw [u^{x, Q2) - d^{x, Q^) 



(27) 



2A 



u^{x,Q^)-d^{x,Q^] 
u^{x,Q^) + d^{x,Q^) 



2x 2 



x,Q' 



(28) 
(29) 



Substituting these into Eqs. (|T^ and (|TH)l and forming the PW ratio in Eq. ([T]), we 

arrive at 



Rj{x,Q' 



2x 



w 



Uy dy 



'2 ^"^'['^ A 3-6xwu^ + d^ 



AA l + (l-y)^<-4 
A l-(l-y)2^^ + 4 



,(30) 



where the arguments x, are implicit in the nPDFs on the r.h.s. 

For isoscalar nuclei, A^ = 0, one obtains the conventional result, R^. = 1/2 — 
Xw For non-isoscalar nuclei, isospin corrections proportional to A^/A arise. We also 
immediately see that in the PW ratio of differential cross sections, i?^(x, Q^), the 
effects of nuclear valence quark modifications cancel out if the modifications are the 
same both for the nuclear valence-w and for valence-d quark distributions. 

We compute the PW ratio in two different ways: 



1. For Rjix,Q 



2. ^NuTeV 



0.2277 and take the same modification Ry from EKS98 for both R^^ and R^^. 



we use the NuTeV value 



X 



NuTeV 
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2. For -R^(x, Q^] {xw), Rdv 7^ ^PP^Y the world average value (xiy) = 0.2227 

and assume R^^ 7^ R^^. 

Setting these two PW ratios equal should then give us an idea into which direction the 
nuclear modifications R^^ and R^^ should deviate from the average valence modifica- 
tion Ry so that the change in the PW ratio caused by the larger NuTeV value of xw 
is accounted for. 

Since the y values we shall consider below, are always larger than 0.1 and Aa/A is 
of the order 0.1 for iron and {uy — dy) / {uy + dy) -C 1, we can expand the denominator 
in Eq. (jSHl)- This gives 

R2{x,Q^,xw) ~ -xh/)|i - (31) 

where h{y, xw) = \'^\[jy^2 — g^g^^ is always positive and not sensitive to small changes 
of Xw, since the first term dominates. Now, we set 

R2ix, Q'; X™, Rl = Rl = R^) = R^^x, Q'; {x^), < ^ <) (32) 

which, since xj^ "^^^ > {xw), leads to 

uy - dy - d^ 

Uy + dy Uy + dy 

where the nuclear effects on the l.h.s. have cancelled each other. Then, applying the 
definitions in Eqs. (j^ - ipUl) to the nuclear PDFs indicates that 

Ri{x,Q') < R^{x,Q') < Rt,{x,Q') (34) 

in the x, region that the NuTeV cross sections are sensitive to. 



4 Results 

Let us now turn to a more detailed quantitative study of the PW ratio Ra{x) in 
Eq. (j2I). In the NuTeV experiment P], the neutrino and antineutrino cross sections 
were measured over a range of beam energies, virtualities and Bjorken-x. Also cuts 
in the calorimeter energy were done. We do not attempt to simulate the kinematics 
in the experimentally collected event sample here in detail. Instead, we simply fix the 
beam energy to -E = 116 GeV, close to the NuTeV's average neutrino and antineutrino 
beam energies and the virtuality to = 20.5 GeV^, corresponding to the NuTeV's 
average virtuality and integrate over x in forming the PW ratio in Eq. To mimic 
the NuTeV cuts on the calorimeter energy, we require the final state energy on the 
hadron-remnant side to be in the range -Emin < m + yE < -Emax, with £^min = 20 GeV 
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and -Emax = 180 GeV. This, together with the requirement y 
the integration hmits for 



Q^/{2mEx) < 1, fixes 



0.094 



2mE 



< X < Xr, 



2m(Emin - m) 



0.57. 



(35) 



We emphasize, however, that exact kinematic hmits are somewhat irrelevant from the 
point of view of this paper. We aim to study whether order-of-magnitude changes in the 
nuclear modifications would be needed in order to explain the NuTeV anomaly in a very 
conventional manner. We shall show that in fact surprisingly modest modifications, of 
the order 20-30 %, are sufficient. In the following, all the results are for iron nucleus 
A = 56 and Aa = 4. 

In Fig.Qwe show the differences of the double differential neutrino and antineutrino 
cross sections from Eqs. (|T^ and (|THjl as a function of computed with (x^/), and 
CTEQ6L [19 PDFs with Rj^ = R^^ = R^ from EKS98 [Hj. From the result, it is 
obvious that the smallest-x region dominates in the x-integrated cross sections. Thus, 
we should be looking for modifications of R^^^ and R^^ in the region x > Xmin ~ 0.1. 
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Figure 1: The cross section differences / dxdQ 



/ dxdQ^ as a function of x computed from Eqs. (|14|) and (|15|) at = 20.5 GeV and 



da^^^ /dxdQ'^ and da^/^'^ / dxdQ'^ 



E = 116 GeV, with EKS98 nuclear PDF effects for iron, A = 56, Z 
world average value of mixing angle {xyi/) = 0.2227. 



26, and using the 



To get transparent order-of-magnitude estimates of the individual modifications 
R^^ and R^^ in light of the NuTeV sin^ 6w anomaly, we assume the following very 
rough form for the valence-d quark modification in iron at = 20.5 GeV^ (see Fig. El 
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ahead) 



{Ci = constant, when x < Xmin 

C2 = constant, when x^m < x < x^^y, (36) 

R^{x, Q2) fromEKS98 when x^ax <x<l 

For a fixed constant C2, the constant Ci is determined from charge conservation 
dxdy{x,Q'^) = 1, after which the ratio R^^ is obtained from i?y(EKS98) and Rj[^ 
according to Eq. (j26|) . Then, as the EKS98- result for Ry conserves baryon number, 
charge conservation QdxUy{x,Q^) = 2 is automatic. We proceed as outhned in the 
previous section. By using the x-integrated cross sections, we compute the PW ratios 
according to Eq. (j21), on one hand for x^y """^^ and setting i?^^ = R^^ = Ry, and on the 
other hand for (xw) and setting R^^ 7^ ^uy- '^^^ requirement (analogous to Eq. (1221)) 

RJiQ^, x^^""^, -Rrf^ = Ruy = Ry) = Ra{Q^^ {^w),Rdy 7^ (2'^) 

now fixes the constant C2 ~ 0.708 in R^^, and consequently Ci ~ 1.16. Figure |21 shows 
the PW ratio as a function of the weak mixing angle, computed with the mutually 
different nuclear modifications Rf and R^ (r.h.s. of Eq. ()37|) . "extreme"), and with 
the EKS98 modification R^ only [l.h.s. of Eq. "EKS98"). 
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Smooth modifications 
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A = 56 ; Z = 26 
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Figure 2: The Paschos-Wolfenstein ratio as a function of the weak mixing angle, computed 
for iron, A = 56, Z = 26, by identifying the individual nuclear modifications R^^ and R^^ 
with the EKS98 Ry, with those in Fig. El ("extreme"), and with those in Fig. El ("smooth"). 
The world average (xw) and the NuteV result indicated. 



2 qNuTgV 



Fig. 121 showing the extreme case where the mutually different nuclear modifications 
to uy and dy explain the whole NuTeV anomaly, is the main result of this paper. Per- 
haps somewhat surprisingly, our result shows that even if we do not invoke any exotic 
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Figure 3: An order of magnitude estimate of the nuclear valence quark modifications needed 
to alone account for the large value of xw reported by NuTeV. The individual modifications 
R^^ and Rf^ are shown for iron, A = 56, Z = 26, as functions of x at a fixed scale = 
20.5 GeV^. The average valence quark modification Ry is from the EKS98-parametrization 
[T^ . Definitions of the ratios Rf are given in Eqs. H23() - l|25() . 
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{x\y) would be explained. 
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new phenomena, such as isospin-breaking PDFs, locally non-zero s — s distributions, 
not to mention beyond-Standard-Model physics, to explain the NuTeV anomaly, fairly 
modest, ~ 20 — 30 % changes in the nuclear valence quark modifications would seem 
to explain the large value of xw observed by NuTeV. 

In addition to the extreme case considered above, we also consider, as an example, 
more modest changes with a continuous and smooth parametrization for Rf^ shown 
in Fig. m Also here charge is conserved and the EKS98 ratio Ry is reproduced as 
explained above. Proceeding again in the same manner as above, we find that these 
nuclear effects would correspond to xw = 0.2243, which would explain a third of the 
NuTeV anomaly Axw = 0.005. 



5 Effects on i?^^ ^dy 

Next, we show that the modifications R^^ and Rf^ obtained above do not induce severe 
modifications into the quantities which are already rather well constrained in the global 
DGLAP analyses of nPDFs. As mentioned above, the data constraints are given by 
the I A DIS cross sections and Drell-Yan dilepton cross sections in pA collisions. 
In DIS of leptons off the free proton, the leading order structure functions are 



F^^{x,Q')=j:elx 



2xF[%x,Q''). (38) 



q{x,Q') + q{x,Q') 
The nuclear modifications of F!^"^ are specified relative to deuterium, 

nU.,Q^)^'^^^^.m^, (39) 



Then, as discussed in ^H], and using the definitions of nPDFs given here, the nuclear 
modifications of F!^^ relative to deuterium become (suppressing the arguments x, Q"^) 



A _ 5(ma + dA + UA + (Ia) + 4(5^ + Sa) + ■■■ + ^(c^A + dA-UA- Ua) 

~ 5(m + d + u + J) + 4(s + s) + . . . 



where the dots denote the heavy quark contributions and the small nuclear effects 
in deuterium have been neglected. From above, we see that since ua + dA = Uy + 
dy + Ua + d-A = Rvi^v + dy) + ua + d-A, any change in R^^ and R^^ that leaves 
the valence modification Ry unchanged, will not at all change the isoscalar part. We 
emphasize that the DIS I A data has usually been approximately "isospin-symmetrized" 
by the experimental collaborations. E.g. the EKS98 analysis is using such isospin- 
symmetrized DIS data (NMC, SLAG mainly) and therefore in constraining the fits, 
one is only discussing isospin symmetric nuclei for DIS. We show in Fig. El the effects 
that our extreme modifications of Rf^ and R^^ in Fig. |21 would induce into the ratio 
Rp^ of non-isoscalar iron at = 20.5 GeV^. Since the changes due to R^^ ^ R^^ only 
take place in the subleading correction term proportional to A a/ A, the effects remain 
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Figure 5: The nuclear modification ratio Rp for non-isospin-symmetrized iron, ^ = 56 and 
Z = 26, as a function of x at a fixed scale Q = 20.5 GeV^. The dashed curve is computed 
with the EKS98 modifications and the solid one is the result when the extreme modifications 
of Fig. I^lare applied for R^^ and R^^- The SLAG data shown have been isospin symmetrized 
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j4 = 56 and Z = 26, as a function of the target-side momentum fraction X2 at a fixed scale 
Q2 = 20.5 GeV^. The dashed curve is computed with the EKS98 modifications and the solid 
one is the result when the extreme modifications of Fig. El are used for R^^ and R^^ . The 
data shown are from E772 j20j . 
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very small, below one per cent level for iron. Given the precision of the existing data, 
such small changes can be easily accommodated in the global DGLAP fits. 

Finally, we perform a similar check for the Drell-Yan nuclear modification ratio 
^dy(^2,<5^) (see [S]), which can be written in the LO as follows 



^daP^/dx2dQ'^ 

= {4[ui{u^ + d^) + ui{u^ + 4)] + [di{d^ + u^) + + 4)] + 4si4 + • • -l/^DY 

+ ^{A[u,{d^ - u^) + - u^)] + [d,{u^ - d^) + d,{4 - 4)]}/iVDY, (41) 

where A^dy = 4[ni(u2 + ^2) + Ui{u2 + ^2)] + [^1(^2 + ^2) + di{u2 + ^2)] + 4siS2 + . • • 
and is the invariant mass of the lepton pair, and where we have used the notation 
q^^^ = q(A){xj,Q'^) with j = 1,2, where X2{xi) is the fractional momentum of the 
colliding parton from the target (projectile). Again the quark combination ua_ + d^ 
appears in the leading, isoscalar, term. For isoscalars, any i?y-conserving changes 
made for R^^ and R^^ will not affect the ratio R^y{x2,Q^) at all. Contrary to the 
DIS case, however, the DY data available are not isospin-symmetrized but show also 
some isospin effects. The modifications of R^^ and Rf^ will be transmitted to R^y 
only through the correction term proportional to A a/ A and hence even the extreme 
modifications of Eq. ()36|) and Fig. El cause only a small deviation from the result 
computed with the EKS98 nuclear modifications. This is shown in Fig. IHl where we 
compute the ratio -Rdy(^25<5^) as a function of X2 at a fixed scale = 20.5 GeV^ 
on one hand with the extreme modifications for R^^ and R^^ and on the other with 
the EKS98 modifications only. Again the difference to the EKS98-based result is less 
than a percent. We thus conclude that the DIS and DY cross sections used in the 
global DGLAP fits are essentially insensitive to the decomposition of Ry into R^^ 
and Rf^ , as long as the average valence quark modification remains unchanged. Thus, 
the interpretation of the NuTeV sin^ 9w anomaly in terms of conventional physics, 
mutually different R^ and R^ , does not lead into a contradiction with the existing 
global DGLAP fits for the uPd'fs. 



6 Effects on R""^ in NOMAD 

As the final point in this paper, we make the following prediction for the NOMAD 
neutrino experiment at CERN ^1]. In this experiment one can measure the ratio of 
charged current and neutral current neutrino-iron cross sections, 

i/,NC 

- %c- (42) 

The differential cross sections are given in Eqs. Q and (|7|). With the parameters 
= 13 GeV^ and {E^) = 45.4 GeV, the typical NOMAD x-range becomes x > 0.15. 
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This is close to that of NuTeV, and the valence quarks dominate the x-integrated cross 
sections. As shown in Fig. [71 we compute the x-integrated cross sections at a fixed 
scale = 13 GeV^ and form the ratio again in two different ways: on one hand 
we apply the EKS98 nuclear effects, where R^^ = R^^ = Ry, and on the other we 
apply the individual modifications R^^ ^ R^^ from Figs. El and El (and EKS98 for the 
parton fiavours other than valence quarks). The figure demonstrates that the situation 
with R\ in NOMAD should be very similar to that with the PW ratio in NuTeV: if 
identical nuclear corrections are used for uy and dy, the value of xw extracted from 
the NOMAD data will be larger than the world average {xw)- If the whole NuTeV 
anomaly is explained by the individual, mutually different, nuclear effects in uy and 
dy distributions, then solving from 

NOMAD p^,EKS98n _ TDUf^2 / ^ \ tdA _ pA,EKS98 / pA n fAo\ 

where the individual modifications R^^ and R^^ are the extreme ones from Fig. |2l 
gives x^^^^ ^ 0.2265, i.e. AxJ^^^^^ = x^^^^ - {xw) ~ 0.004, close to the result 
obtained by NuTeV. With the smooth modifications in Fig. |31 the increase in xw is 
again close to that obtained in Fig. |21 
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Figure 7: The ratio i?^ for iron, A = 56, Z = 26, as a function of sin^ 9]y, computed in the 
x-range of NOMAD at = 13 QeV^ with the EKS98 nuclear modifications ("EKS98"), with 
extreme individual nuclear modifications R^^ and R^^ in Fig. O together with the EKS98 
for other parton flavours (" extreme" ) , and with the smooth modifications in Fig. jl] together 
with the EKS98 ("smooth"). The world average (xw) and the expected increase of xw iu 
the NOMAD case, cf. Eq. (HHJ, are indicated. 
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7 Conclusions 



For the nuclear valence quark PDFs one usually makes the approximation Uy/uy ~ 
dy/dy ~ {uy + dy)/{uy + dy) = Ry. In this case, the effects from the nuclear valence 
quark modification Ry are cancelled out in the Paschos-Wolfenstein ratio R~{x,Q^). 
In this paper, we are suggesting that taking R^^ = Uy/uy ^ dy/dy = R^^ and by 
allowing nuclear effects of the order of 30 %, leads to the modifications of the Paschos- 
Wolfenstein ratio R'iQ"^) which are as large as would be induced by increasing the 
weak mixing angle xw to the large value reported by NuTeV. We have also checked 
that as long as the valence modification Ry is conserved in the decomposition into 
R^^ and -R^^, the computed ratios of the nuclear DIS and DY cross sections over 
deuterium, used as data constraints in the global DGLAP analyses of nPDFs, change 
by less than a percent in the case of iron. On one hand, this insensitivity shows that the 
explanation of the NuteV anomaly in terms of nPDFs is well possible without running 
into contradiction with the global nPDF analyses. On the other hand, it shows, in 
agreement with [131111111, that in practice it is not possible to pin down the individual 
nuclear effects R^^ and R^^ based on finite-precision nuclear DIS and DY data. The 
main point of the present paper thus is that most, if not all, of the siv? 6w NuTeV 
anomaly could be explained by conventional physics and that the NuTeV result can 
play a key role in disentangling these effects for uy and dy. Finally, as a consequence, 
we predict that if the whole NuTeV anomaly can be accounted for by the nuclear 
corrections as we suggest here, and if identical nuclear effects for uy and dy are used 
in the data analysis, the value of sin^ 6w to be determined in the NOMAD neutrino 
experiment at CERN should be quite close to that obtained by NuTeV. 

We wish to emphasize that in the global nPDF analyses the shapes of the individual 
nuclear corrections to the PDFs are determined based on the data constraints (through 
DGLAP evolution) and sum rules. The origin of these corrections, and PDFs in general, 
is nonperturbative. If the NuTeV anomaly indeed is due to the nPDF effects as we 
suggest above, it would be very interesting to understand the mechanism which makes 
the nuclear effects for uy and dy distributions so different. We do not have an answer 
to this question for the moment. 

In this paper, we do not make an attempt to simulate the NuTeV kinematics in 
detail but our goal is to show that interpretation of the NuTeV anomaly in terms of 
nPDF effects is indeed possible and to see the order of magnitude for the effects needed. 
The more detailed kinematics, along with scale evolution details of the valence quark 
effects, we shall consider elsewhere. 
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